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1.   Introduction. 

Motivated  by  applications  to  reliability  theory,  the  model  we  shall 
consider  is  a  mixture  of  the  form 

oo 

FT(t)  =    e"XtdG(X),     t  £  0  (1.1) 

0 

where  F  (t)  is  the  survival  function  of  the  mixture,  F(t;A)  =  e  is  the 
survival  function  of  the  component  distribution  given  that  A  =  A  for  A  a 
nonnegative  random  variable  having  cdf  (cumulative  distribution  function)  G. 
We  can  refer  to  such  a  mixture  as  a  G-mixture  of  exponentials.  In  the  relia- 
bility context,  T  represents  the  time  to  failure  of  an  item  drawn  at  random 
from  a  population  of  items  with  different  exponential  life  distributions,  mixed 
according  to  the  distribution  of  A. 

We  take  the  usual  definition  for  a  failure  rate  function,  i.e.  for  the 
random  variable  T  having  density  f„(t)  and  survival  function  F  (t) ,  the 
failure  rate  function  is  defined  by 

fT(t) 

h  (t)  =  - ,    for  all  t  such  that  F  (t)  >  0. 

FT(t) 

The  initial  mixing  distribution  G  may  be  modified  to  yield  a  related 
distribution  at  any  time  t,  i.e.  let  G  (A)  =  P(As.A  |  T>t)  and  define  A(t) 
to  be  a  random  variable  having  cdf  G  .  G  is  the  residual  mixing  distribu- 
tion at  time  t.  It  represents  the  revised  distribution  for  the  failure  rate 
of  a  randomly  drawn  item  from  the  mixture,  given  that  the  item  is  still  alive 
at  time  t.  The  relationship  between  G  and  G  may  be  exhibited  through  a 
straightforward  application  of  Bayes'  Theorem  which  yields 


-At 

dG.(X) ^r—  dG(X),  (1.2) 

Z  E(e  At) 

where  the  expectation  is  with  respect  to  G,   the  distribution  of  A. 

For  such  a  mixture  of  exponentials  as  (1.1),  it  is  well  known  (Barlow, 
Marshall,  and  Proschan  [2])  that  the  mixture  failure  rate  h  (t)   is  decreasing  in 
(we  shall  use  'decreasing'  for  non-increasing  and  'increasing'  for  non-decreas- 
ing), and  that  (Gnedenko,  Belyayev,  and  Solovyev  [3])  the  mixture  failure  rate 
function  approaches  as  a  limit  with  increasing  time   t   the  least  of  all 
parameter  values  positively  present  in  the  mixture,  i.e.  X    . 

Aldrich  and  Morton  [1]  and  O'Bar  [4]  showed  that  at  time  t   the  failure 
rate  function  of  a  mixture  of  exponentials  is  the  expected  value  of  the  resid- 
ual mixing  distribution  corresponding  to  that  time,  i.e. 

hT(t)  =  E[A(t)]  (1.3) 

Thus,  the  sequence  of  residual  mixing  distributions  is  decreasing  in  expectation. 
This  fact  supports  the  intuitive  notion  that  with  increasing  time  the  shorter- 
lived  members  of  the  mixture  are  being  "weeded  out." 

The  conjecture  of  a  compatible  stochastic  ordering  of  the  residual 
mixing  distributions  arises  as  a  logical  extension  of  the  foregoing  results. 
In  this  paper  the  existence  of  such  a  stochastic  ordering  is  demonstrated.   It 
is  also  shown  that  the  residual  mixing  distributions  converge  in  distribution 
to  the  distribution  degenerate  at  the  least  parameter  value  positively  present 
in  the  mixture,  with  an  attendant  monotone  convergence  of  raw  moments.   The 
latter  result  is  one  asserted  by  Aldrich  and  Morton. 


The  results  are  presented  as  applications  of  Propositions  1  and  2,  to 
be  developed  in  the  next  section.   These  propositions  state  properties  of 
E{u[A(t)]}   for  fairly  general  functions  u.   To  indicate  the  relevance  of 
these  propositions,  we  first  note  that  an  expression  of  this  form  is  involved 
in  (1.3),  taking  u[A(t)]  =  A(t).   We  next  define  an  indicator  function 


I,  (A) 


J 


o, 


1, 


A^A 


A  >  A 


(1.4) 


Using  this  indicator  function,   G   may  be  expressed  as 


GJX)  = 


dGt(x)  =  E{Ix[A(t)]}, 


(1.5) 


another  expression  of  the  same  form,  with  u[A(t)]  =  I.[A(t)] 


2.   Basic  Results . 

We  shall  develop  in  Propositions  1  and  2  of  this  section  basic  tools 
which  should  prove  to  be  useful  in  a  wide  variety  of  applications.  Several 
of  these  applications  are  given  in  section  3. 

Using  (1.2)  we  find  that 


E{u[A(t)]}  = 


u(A)dGt(A) 


'  u(A)e 


-At 


E(e-At) 


dG(A) 


m   E[u(A)e"At] 
E(e-At) 


(2.1) 


Thus,  expectations  of  functions  of  A(t)   may  be  evaluated  by  simply  taking 
expectations  of  functions  of  A.   A  straightforward  application  of  the  Monotone 
Convergence  Theorem  yields 

£  {E[u(A)e"At]>  =  E[u(A)  ±   e_At] 

=  -E[Au(A)e"At],  (2.2) 

for  A  £  0. 

We  next  require  the  function  u  to  be  nonnegative  and  monotone  to  get 

PROPOSITION  1.   If  A  is  a  nonnegative  random  variable  and  u  a  nonnegative 
monotone  increasing  (decreasing)  function,  then 

.  E[u(A)e-At]  .  E{u[A(t)]} 

E(e  At) 

is  monotone  decreasing  (increasing)  in  t. 

Proof:   Using  (2.2), 

A.  -  (t)   m   E(e"At)E[-Au(A)e"At]  -  E[u(A)e"At]Er-Ae"At] 
dt  V  E(e-At)2 

E[Au(A)e"At]    E[u(A)e"AT]   E(Ae"At) 
E(e"At)        E(e-At)   '  E(e"At) 

=  E{u[A(t)]}E[A(t)]  -  E{A(t)u[A(t)]} 

=  -  C0V{u[A(t)],  A(t)}. 

Then  by  the  proporty  of  the  covariance  of  similarly/dissimilarly  ordered 
functions  of  the  same  random  variable, 


(i)   u  increasing  in  A(t)  =>  —  (p  (t)  £  0,   and 
(ii)   u  decreasing  in  A(t)  »  —  <p    (t)  ^0.  Q 

We  now  need  two  lemmas  regarding  limits  of  expectations  involving 
time   t  as  a  parameter. 

LEMMA  1.   Given  u  a  nonnegative  function, 

lim  E[Ix(A)u(A)e"(A_X)t]  =  0. 

Proof:   Note  that  I.(A)=1*»A-X>0*»e"^+  in  t.   An  application 
of  the  Monotone  Convergence  Theorem  yields  • 

lim  E[Ix(A)u(A)e"(A_X)t]  =  E[I  (A)u(A)  lim  e~(A~X)t]  =0.      D 
t-x»  t-*a> 

LEMMA  2.   For  X  >  X  =   inf{X  I  G(X)<1}, 
o        ' 

lim  E[(l-I.  (A))e"(A"X)t]  >  0. 

A 
t-x» 

Proof:      Note    that      (1-1    (A))    =1*>A-A*0*>  e-(A_x)t   fe  ±      for      fc   B  Q;      Then 

A 

E[(l-I .  (A))e"(A_X)t]  I>E[1-I.  (A)]  =  G(X),   and   lim  E  [  (1-1,  (A)  )e~(A_X) C ]  £  G(X)  >  0 
A  A  A 

t-x» 

D 

If,  in  addition  to  the  mono tonicity  requirement  of  Proposition  1,  we 
also  require  u  to  be  right-continuous  at   X  ,   we  may  find  the  limiting  value 
of  cp  (t)   as   t  •*■  •.   We  state  this  result  in 

PROPOSITION  2,   Let   A  be  a  nonnegative  random  variable  having  cdf   G,   and 
u  a  nonnegative  monotone  increasing  (decreasing)  function.   If  u  is  right-  - 


continuous  at  X  =  inf{X  |  G(X)<1},   then  <p  (t)   decreases  (increases) 

monotonically  to  u(X  )   as  t  -*■  °». 

o 

Proof: 

(i)   If  u(A)t   in  A  then  cp  (t)  ^  u(X  ).   By  Proposition  1,  <p  (t)+  in  t. 

Thus,   lim  u>  (t)   exists  and  lies  in  the  interval   [u(X  )  ,°°)  .   Choose 
u  o 

X  >  X  .   Then, 
o 

E[(l-Ix(A))u(A)e"At]  +  E[Ix(A)u(A)e"At] 

(p  (t)  =  jr ^7 

U         E[(l-lA(A))e  At]  +  E[Ix(A)e  At] 

E[(l-Ix(A))u(A)e"(A"X)t]  +  E[Ix(A)u(A)e"(A_X)t] 


E[(l-Ix(A))e_(A-X)t:]  +  E[Ix(A)e-(A_X)t] 


Since  u(A)+   in  A 
cp  (t)  £> 


u(X)E[(l-Ix(A))e  (A  X)t]  +  E[Ix(A)u(A)e  (A  X)t] 


[Ix(A)e 
-(A-X)t 


E[(l-Ix(A))e  v"  """]  +  E[I,(A)e  v"  "/l-] 


E[Ix(A)u(A)e  v"  "'"] 
E[(l-Ix(A))e  (A  X)t]  +E[Ix(A)e  (A  X)t] 

Lemmas  1  and  2  may  be  applied  to  show  that  the  last  term  in  expression 

(1)  goes  to  zero  as   t  -*■  °°.   Therefore,   lim  <p  (t)  ^  u(X)  .   Let  X  4-  X 

to  get   lim  cp  (t)  £  u(X  ).   Then,  cp  (t)  +  u(X  )   as   t  ■>  °°. 
°         u         o  u         o 

t-x» 

(ii)   The  proof  for  u  monotone  decreasing  is  similar  to  that  above.    D 


3.   Applications. 

The  first  application  of  the  tools  developed  in  the  preceding  section 
concerns  the  stochastic  ordering  of  the  sequence  of  residual  mixing  distribu- 
tions and  its  limiting  distribution,  stated  as 


Theorem  1.   If  A  is  a  nonnegative  random  variable  with  cdf  G  and 

X  =  inf{X  I  G(X)<1},   then  for  a  G-mixture  of  exponentials 
o        ' 

a)  G   is  a  stochastically  decreasing  sequence  in  t,   i.e. 

t  &  t  »  G   (X)  £  G   (X),   for  all  X. 
1        2 

b)  G   converges  in  distribution  to  the  distribution  degenerate  at 

X  ,   as  t  ■>  °°. 
o 

Proof: 

a)  By  (1.5), 

Gfc(X)  =  E{lA[A(t)]>. 

Note  that  I,   is  monotone  increasing  in  its  argument.   Application 
of  Proposition  1  with  u[A(t)]  =  I,  [A(t)]   yields  G  (X)4-   in  t. 

A  t 

b)  Application  of  Proposition  2  with  u[A(t)]  =  I, [A(t)]   yields  as 

t  ■>  °°, 


0,   X  £  X 
o 


Gt(X 


(  l,  x  >  x  .  n 


An  immediate  consequence  of  Proposition  2  is  the  well-known  result  that 
a  mixture  of  exponentials  is  DFR,  with  its  failure  rate  approaching  the  least 
parameter  value  positively  present  in  the  mixture.   We  state  this  result  as 

Theorem  2.   If  A  is  a  nonnegative  random  variable  with  cdf  G  and 

X   =  inf{X  |  G(X)<1},   then  for  a  G-mixture  of  exponentials  the  mixture  failure 

rate  function,  h  (t),   decreases  monotonically  to   X   as   t  -*■  °°. 


Proof:   As  noted  in  (1.3), 


hT(t)  =  E[A(t)]. 


Application  of  Proposition  2  with  u[A(t)]  =  A(t)   yields  the  desired  result. 

D 
We  may  also  see  the  monotone  convergence  for  the  raw  moments  of  the 
residual  mixing  distribution  in 

Theorem  3.   For  p  £  0,   E[A(t)  ]   decreases  monotonically  to   A  . 

Proof:   Application  of  Proposition  2  with  u[A(t)]  =  A(t)   yields  the  desired 
result.  □ 

The  variance  of  A(t)   is  of  interest,  since  it  has  an  interpretation 
as  the  negative  of  the  slope  of  the  mixture  failure  rate  function  (O'Bar  [4]). 
We  state  the  following  as 

Corollary  1.   The  limit  as   t  ->  °°  of  VAR[A(t)]   is  zero. 

Proof:   We  express  the  variance  of  A(t)   as 

VAR[A(t)]  =  E[A(t)2]  -  E2[A(t)]  (3.1) 

Then,  by  Theorem  3, 

lim  VAR[A(t)]  =  A2  -  A2  =  0.  □ 

o    o 

t-Hn 

Expression  (3.1)  leads  us  to 

Corollary  2.   On  the  interval   (0,°°),   VAR[A(t)]   is  of  bounded  variation. 

Proof:   From  the  power  series  expansion  for  e   ,  we  have  that 

An   n!   At 
A  <  —  e   , 
n 
t 


for  n  a  nonnegative  integer.   Thus, 


„,.n  -At,     ,  „r  At  -At, 
E[A(t)n]  =  E'Ae„  '  <  4  Ele  _*„   1 


E[e-At) 


n  „/  -Atx 
t   E(e   ) 


n! 


nw  -Atv 
t  E(e   ) 


,   for  t  >  0. 


Then,  from  expression  (3.1)  and  Theorem  3,  since  VAR[A(t)]  is  expressed  as 
the  difference  of  two  real-valued  nonnegative  monotone  functions  for  t  >  0, 
it  is  of  bounded  variation  on   (O,00).  D 
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